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Abstract 

Many of the hadron-hadron cross sections required for the study of the dy- 
namics of matter produced in relativistic heavy-ion colhsions can be calculated 
using the quark-interchange model. Here we evaluate the low-energy disso- 
ciation cross sections of J/ip, tp' , Xj '^^ ^^d T' in collision with vr, p, and 
K, which are important for the interpretation of heavy-quarkonium suppres- 
sion as a signature for the quark gluon plasma. These comover dissociation 
processes also contribute to heavy-quarkonium suppression, and must be un- 
derstood and incorporated in simulations of heavy-ion collisions before QGP 
formation can be established through this signature. 



I. INTRODUCTION 



The first collisions of heavy-ion beams at the Relativistic Heavy-Ion Collider (RHIC) at 
Brookhaven National Laboratory heralded a new era in the study of matter in the extreme 
conditions of very high temperatures and energy densities. The ultimate goal of this research 
is the production of a quark-gluon plasma, which is an unusual phase of strongly interacting 
matter that purportedly existed shortly after the Big Bang [^0. 

The search for the quark-gluon plasma relies on the unusual properties of the plasma 
for its detection. However, many conventional hadrons are also produced during a heavy 
ion collisions. Whatever signal is chosen for the identification of the quark-gluon plasma, 
contributions to that signal from conventional hadronic processes must be identified as back- 
grounds and removed from the data. 

Although there are recent, tantalizing hints of possible quark-gluon plasma production 
in heavy-ion collisions at CERN ||^,^, conclusive evidence is still lacking due to uncertain 
backgrounds from conventional hadronic sources. Investigations of the various 'hadronic 
background processes' are urgently needed if we are to develop a satisfactory understanding 
of the various signals proposed as signatures of a quark-gluon plasma. In this paper we 
consider one type of hadronic background processes that can contribute to heavy- quarkonium 
suppression, which is frequently cited as a QGP signature. 

Matsui and Satz originally suggested the use of suppressed J/ip production as a 
signature for the formation of a quark-gluon plasma in high-energy heavy-ion collisions. 
The recent experimental observation of anomalous J / ip suppression in Pb-|-Pb collisions by 
the NA50 Collaboration has been considered by many authors |]8Hl7|. However there is 
considerable uncertainty in these studies, due to the lack of reliable experimental information 
on J / ip and xj dissociation cross sections in low-energy collisions with light hadrons. Because 
heavy quarkonia decay strongly, many of the dissociation cross sections cannot be measured 
directly in hadron scattering experiments; the cross sections are instead typically estimated 
using theoretical models. Evaluation of these cross sections is of particular interest for 
clarifying the physics of the J/ip anomalies observed in Pb-Pb collisions, and may be of 
considerable importance in future J /ip studies using the RHIC and LHC colliders. 

The dissociation of the J/ip by hadrons has been considered previously in several theo- 
retical studies, but the predicted cross sections show great variation at low energies, largely 
due to different assumptions regarding the dominant scattering mechanism |lT8|- p^] . 

Kharzeev, Satz, and collaborators [^,^ employed the parton model and perturbative 
QCD "short-distance" approach of Bhanot and Peskin [p9| , p0| , and found remarkably small 
low-energy cross sections for collisions of J / with light hadrons. For example, their J / ip+N 
cross section at ^/s = 5 GeV is only about 0.25 [|18|. A finite-mass correction increases 



this cross section by about a factor of two [0 . However, in high-energy heavy-ion reactions 



the collisions between the produced vr and p with J/ ip and ip' occur at low energies (typically 
from a few hundred MeV to about 1 GeV relative kinetic energies) . The applicability of the 
parton model and pQCD for reactions in this low-energy region is certainly open to question. 
Matinyan and Miiller Haglin p5[, Lin and Ko [Q, and Oh, Song, and Lee re- 



cently reported results for these dissociation cross sections in meson exchange models. These 
references all use effective meson Lagrangians, but differ in the interaction terms included 
in the Lagrangian. Matinyan and Miiller included only t-channel D meson exchange, and 



2 



found that the dissociation cross sections of J/ip by vr and p are rather small; both are 
~ 0.2-0.3 mb at = 4 GeV. Including form factors (arbitrarily chosen to be Gaussians 
with a width set to 1.5 GeV) would reduce these cross sections by an order of magnitude. 
Haglin obtained a very different result, with much larger cross sections, by treating the 
D and D* mesons as non-Abelian gauge bosons in a minimally coupled Yang-Mills meson 
Lagrangian. Form factors were also introduced in these calculations ||25|-^. The resulting 
mb-scale cross sections are very sensitive to the choice of form factors. Charmonium dis- 
sociation by nucleons has also been considered recently using a similar effective Lagrangian 
formulation [^. Of course the use of a Yang-Mills Lagrangian for charmed mesons has 
no a priori justification, so this crucial initial assumption made in these references requires 
independent confirmation. The assumption of the t-channel exchange of a heavy meson such 
as a D or D* between a hadron and a J/ip is also difficult to justify physically, because the 
range of these exchanges (1/M ^ 0.1 fm) is much smaller than the physical extent of the 
initial hadron and the J/ip. 

Charmonium dissociation processes can undoubtedly be described in terms of interquark 
interactions, as we attempt in this paper. Since these reactions are of greatest phenomeno- 
logical interest at energy scales in the resonance region, we advocate the use of the known 
quark forces to obtain the underlying scattering amplitudes from explicit nonrelativistic 
quark model hadron wavefunctions of the initial and final mesons. 

Martins, Blaschke, and Quack previously reported dissociation cross section calcu- 
lations using essentially the same approach we describe here. The short-distance interaction 
used by these authors in particular is quite similar to the form we employ. For the confining 
interaction, however, they used a simplified color-independent Gaussian potential between 
quark- ant iquark pairs, rather than the now well-established linear X{i) ■ X{j) form. They 
found a rather large Tc+J/ip dissociation cross section, which reached a maximum of about 
7 mb at a center-of-mass kinetic energy Eke of about 0.85 GeV. Although our approach 
is very similar to that of Martins et ai, our final numerical results differ significantly, due 
mainly to our different treatments of the confining interaction. 

In this paper we use the approach discussed above to evaluate the dissociation cross 
sections of J/ip, ip', Xi cind T by vr, p, and and compare our results to other theoretical 
cross sections reported in the literature. The dissociation cross sections of xj mesons are of 
special interest, as about 1/3 of the J/ip mesons produced in a high-energy nucleon-nucleon 
collision come from the decay of x states 0]. The dissociation cross sections for T are 



also interesting and they have recently been estimated to be quite small in an effective- 
Lagrangian meson exchange model because of large thresholds for the dissociation of T by 
both TT and p 

We employ the Barnes-Swanson quark-interchange model [p3| , p4| to evaluate these dis- 
sociation cross sections in terms of wavefunctions and interactions at the quark level. We 
use the nonrelativistic quark potential model and its interquark Hamiltonian to describe the 
underlying quark-gluon forces. The model parameters and quark masses are determined by 
the meson spectrum, so there is little additional freedom in determining scattering ampli- 
tudes and cross sections. We thus imphcitly incorporate the successes of the quark model 
in describing the meson spectrum and many static and dynamical properties of hadrons. 
We proceed by calculating the scattering amplitude for a given process at Born order in the 
interquark Hamiltonian; the good agreement of this approach with experimental scattering 
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data on many low-energy reactions [p3|-|36[] provides strong motivation for the application 
of this approach to hadron reactions in relativistic heavy-ion collisions. A brief summary of 
the present work has been reported previously pl[] . 

This paper is organized as follows. In Section II we summarize the Barnes- Swanson 
model of quark interchange as applied to the calculations of the dissociation cross sections. 
The reaction matrix can be described in terms of the "prior" or "post" diagrams, which 
are discussed in Section III. Section IV gives some details of the evaluation of the spin and 
spatial matrix elements for the general meson-meson scattering problem. In Section V the 
spin matrix elements are derived explicitly in terms of 9-j symbols. The evaluation of spa- 
tial overlap integrals for the case of all S'-wave mesons is discussed in Section VI. In Section 
VII, we evaluate the corresponding overlap integrals for one P-wave meson. An accurate 
determination of these matrix elements requires correspondingly accurate bound state wave- 
functions; the evaluation of these wavefunctions is discussed in Section VIII. The numerical 
agreement between the post and prior scattering formalisms is demonstrated explicitly in 
Section IX, which provides a very nontrivial check of the accuracy and internal consistency 
of our calculations. Section X presents our results for the dissociation cross section of J/ip 
and in collision with various light mesons, and Section XI gives the corresponding cross 
sections for the dissociation of T and T'. Section XII present results for the dissociation 
of P-wave charmonium states, the xj mesons, in collision with tt, p, and K. Finally, we 
present conclusions in Section XIII. 



II. THE MODEL 

We shall briefly summarize the model of Barnes and Swanson for constituent-interchange 
processes in the reaction 

A(12) +P(34) ^ C(14) +P)(32) (1) 

where A, B, C, and D are qq mesons, and 1, 3, and 2, 4 label the quark and antiquark 
constituents respectively. In this meson-meson scattering problem the scattering amplitude 
in the "prior formalism" is the sum of the four quark-line diagrams of Fig. 1. These are 
evaluated as overlap integrals of quark model wavefunctions using the "Feynman rules" given 
in App. C of Ref. |3^. This method has previously been applied successfully to the closely 



related no-annihilation scattering channels 1 = 2 tttt I = 3/2 Kir / = {0,1} 



5'-wave KN scattering and the short-range repulsive A^A^ interaction |37|. 

The interaction Hij{rij) between the pair of constituents i and j is represented by the 
curly line in Fig. 1 and is taken to be 

' {^olor— Coulomb ('"jj) ~l~ linear (^jj) ~l~ ^pin— spin(^ij) ~l~ ^on} 

= M . M K _ 36 _ _ / ^\ .^.^ I 

2 2 [rij 4 ^ SmiTTij ^ \tt^I^ j J' ^ ' 

where is the strong coupling constant, h is the string tension, and irtj are the masses of 
the interacting constituents, a is the range parameter in the hyperfine spin-spin interaction. 
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and V^on is a constant. For an antiquark, the generator A/2 is replaced by — 





C2 





Tl T2 

Fig. 1. 'Prior' diagrams for Born-order meson-meson scattering. 

It is convenient to introduce Vijirij) to denote the quantity in curly brackets in Eq. (H) 
so that we can write Hij{rij) in the form 

^..M = ^-^ V;,M. (3) 

The Born-order T-matrix element Ty, is proportional to the matrix elements hji of this 
residual interaction (as defined in Ref. |3^). For each of the scattering diagrams of Fig. 1, 
hij and Tij are given as the product of four factors, 

hfi (27r)'^ '^■f^ ^ -^flavor -^color -^spin— space • (4) 

The overall sign S" is a fermion-permutation phase known as the "signature" of the diagram, 
which is equal to (—1)^"=, where is the number of fermion line crossings. {S = —1 for 
the diagrams in Fig. 1.) The flavor matrix element /flavor is the overlap of the initial and 
final flavor wavefunctions. In all the processes considered in this paper, /flavor is equal to 1 
for all diagrams. The color matrix element /color is —4/9 for diagrams CI and C2 of Fig. 1, 
and is +4/9 for diagrams Tl and T2. The spatial and spin matrix element /spin-space is the 
matrix element of Vij, and can in general be written as a sum of products of a spin matrix 
element /spin times a spatial matrix elements /space- The spin matrix element /spin involves 
Clebsch- Gordon coefficients and the spins of the colliding particles and is tabulated for all 



cases of 5- wave mesons in |33|. An explicit closed- form expression for this /spin in terms of 



Wigner's 9j symbols will be given in Section V. The evaluation of the spatial matrix element 
/space will be discussed in detail in Sections VI and VII. 

For the reaction A + B ^ C + D, with an invariant momentum transfer t 

t = [A - Cf = m\ + ml - 2AqCo + 2A-C, (5) 

the differential cross section is given by 

MM' (6) 



d<^fi _ 1 I » ^ |2 



12 



dt 647rs|pyi| 
where the matrix element Ai ji is related to Tfi by 
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Mfi = ^{2Ea){2Eb){2Ec){2Ed) Tj^. 



(7) 



In Eqs. (6) and (7), and Ea are the momentum and the energy of meson A in the 
center-of-mass system. The total cross section for the reaction A + B ^ C + D can be 
obtained from dafi/dt by integrating over t. 



III. POST AND PRIOR DESCRIPTIONS 

Before proceeding to our results, we note that a well-known "post-prior" ambiguity arises 
in the calculation of bound state scattering amplitudes involving rearrangement collisions 
||39|| . Since the Hamiltonian which describes the scattering process AB —>■ CD can be 
separated into an unperturbed Hamiltonian and a residual interaction in two ways, H = 



h'^a + Hq' + Vab = Hq' + H^^' + VcD, there is an ambiguity in the choice of Vab or Vcd as 
the residual interaction. The first version gives the "prior" diagrams of Fig. 1, in which the 
interaction occurs before constituent interchange. The second choice is the "post" formalism 
in which the interaction occurs after constituent interchange, as in the diagrams of Fig. 2. 



r(0) 



r(0) 



r(0) 



A 
B 




CI 




C2 





Tl T2 

Fig. 2. 'Post' diagrams for Born-order meson-meson scattering. 

One may prove in the context of non-relativistic quantum mechanics that the 'prior' 
and 'post' diagrams give the same scattering amplitude and hence the same cross section, 
provided that exact bound state wavefunctions of the various are used for the external 

states (This is discussed in detail and is demonstrated numerically in Ref. |Q for 

7171 pp scattering.) A consistent calculation thus leads to description-identical results 
for the scattering amplitude in non-relativistic quantum mechanics. We shall confirm the 
"prior-post" equivalence numerically in our non-relativistic calculations of the J/iJj and ijj' 
dissociation cross sections. 



IV. EVALUATION OF THE MATRIX ELEMENT /spin-space 

For the processes considered here, it suffices to treat reactions of the form ^4(12) + 
5(34) — > C(14) + -D(32), in which constituents (antiquarks) 2 and 4 are interchanged, as 
depicted in Figs. 1 and 2. We denote the total angular momentum, the orbital angular 
momentum, and the spin of meson a {a = A, B, C, and D) by J^, La, and Sa respectively, 
with the associated spatial wavefunction $a and spin wavefunction Xa- 
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The quantity /spin-space is the matrix element of Vij{rij) between the initial and final 
meson states. The interaction Vijirij) is the spin and spatial part of Hij (Eqs. @ and (H)), 
and it consists of the sum 



(8) 



where the superscripts (z) represent the color- Coulomb, linear, and spin-spin interactions 
respectively. Specifically, v^p = t>p) = 1, f^^^ = Sj ■ Sj, and the corresponding t>W can be 
determined from Eqs. and (^). For the scattering problem the sum of the amplitudes of 
all diagrams arising from the constant term V^on is zero, so we do not need to include Kon 
in deriving scattering amplitudes and matrix elements. 



The matrix element /spin-space is therefore the sum of three terms, each of which is of the 



form 



{{^CXc)ZS^DXD)Z>rVs\{^AXA)ZS^BXB)TJ 

= X! {Jc Jcz Jd Jdz\J' J'z){.J Jz\J aJaz Jb Jbz 
JJ^J'J' 



{ [($CXc)-'^(<fDXD) 



Jd 



J' 



i<^AXByH<i>BXBy^] , ). (9) 



In the above equation, the coupled initial state | aXb)'^"^ bXb)'^'^ j ) = l^m ^) of mesons 
A{12) and 5(34) can be written as g^] 

\^-^^')=T.iiXAXBf{^A^ByjJzmAXByH^BXBy^JJz)\ UxAXBy {^A<^By]' ) 

Sa Sb S 
^A Lb 

J A Jb J 



S,L 



^ SUaJb {La Lb L \ {S L (J, - S,)\J J,) | (XAXB)l($A$B)fj,_5.) ), (10) 



SLS^ 



where S = ^2^ + 1. The final state | [($cXc)^^($dXd)^°1 ) = l^^nlf) of mesons C(14) 

and D{32) can be written in a similar way, so the matrix element on the righthand side of 
Eq. (D is 





' Sa 


Sb 






' Sc 


Sd 


E SLJaJb < 


La 


Lb 


;l 


> S'L'JcJd < 


Lc 


Ld 


SLSzS'L'Si 


J A 


Jb 






. Jc 


Jd 



x{SS,L (J, - S,)\J J,){S' Si L' (J^ - 5^)1 J' 



X {{<^c'^Dr(J'^-S'^)\Vr\{^A^Br(J.-S.)) ((XcXD)IJt^. | • (H) 

The above result shows that /spin-space is in general a sum of products of a 
spatial matrix element {{(^c^d)[ji -s')\''^r\{.^A^B)[,j^^s^)) ^ ^Pi^ matrix element 



-^spin 



{{XcXD)s'J^s\iXAXB)sJ- For our interaction, the spin matrix element 
{iXcXD)s[\^s\iXAXB)sJ is diagonal in S and 5*2, and is independent of 5*2, as shown in 
the next section. 
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In this paper, we shall specialize to the cases in which mesons B, C and D are all S'-wave 
mesons with Lb = Lq = = 0. Therefore we have L = La, Jb = Sb, J' = S, J'^ = S^, 
and 

Sa Sb S 

= SLaJaJb {La La } (S S, La (J. - S,)\J J,) 
Ja Sb J 



final I 



{^C^D\Vr\{^A^B)lt^sJ{iXcXD)sJVs\{XAXB)sJ, 



(12) 



where \Sa — Sb\ < •S' < {Sa + 5*^) and \S — La\ < J < {S + Lj 



For the collision of 



unpolarized mesons, we can calculate the square of the matrix element, | /space-spin | , average 
it over the initial states, and sum it over the final states. The result is 



^ space— spin | 



Sa Sb S 
La La 
Ja Sb J 



x\{S S^La {J,-S,)\J J, 



E(*C$D I I {^A^B)'(i-S.)) ( {XcXD)l\vf I {XaXb)1 ) 



■ (13) 



The summation over 5"^ can be carried out and the summation over J^ can be converted to 
a summation over Ma- We then obtain 

s V 



'space— spin | 



Sa Sb 

J2 (SJf {La La 

Ja Sb J 

3 



SJMa 



Y.^^C^D\v'^\^A^B)lt){{XcXD)l;MP\{XAXB)l) 



(14) 



a{LAMASS, 



(15) 



where —La < Ma < La- From the relation between the matrix element of Vij{rij) and 
the cross section, the above result leads to the following "unpolarized" cross section for the 
collision of unpolarized mesons, 

' Q O O ^ 

OA OB o 

La La 
Ja Sb J 

where cILaMaSSz) is the cross section for the initial meson system to have a total internal 
orbital angular momentum La and total spin S with azimuthal components Ma and 5*2 
respectively. For our interaction of Eq. (|^), ^(L^M^S'S'z) is independent of S^, and thus the 
label Sz can be omitted. We can write out the results for other simple unpolarized cases. If 
La 7^ and Sb = 0, then S = Sa and the result of Eq. (p!5| ) becomes 



^unpol^ E {SJf< 
SJMa 



a 



unpol 



1 



La 



(2L, 



1^ 



J2 ct^LaMaSa)- 



Ma=-La 

If = and 5*^ 7^ 0, then the result of Eq. (^) is 

where cr(S') is the cross section when the initial two-meson system has a total spin S. 



a 



unpol 



(16) 



(17) 
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V. EVALUATION OF THE SPIN MATRIX ELEMENT 



We denote the spins of the constituents in the scattering process y4(12)i?(34) — »• 
C(14)D(23) by Si, S2, S3, and S4. Using properties of the Wigner {9j} symbols |^^, we 
may rearrange the spins to obtain 

\{xaXb)sJ = I [(siS2)5a(s3S4)5b]^J 

f Si S2 Sa ] 

= (_1)5b-S4-S3 ^ SaSbSiAIs, S3 5b |[(SiS4)Si4(s2S3)^23]l) (18) 

•S'i4 523 [ 5'i4 ^23 S J 

The phase factor (— i)'5s-«4-s3 arises from an interchange of spins in the Clebsch-Gordon 
coefficients, 

\{s3S,)Sb) = hlf^-''-''\{s,ss)SB). (19) 
The matrix element of the spin unit operator f ^ = 1 is then given by 

{iXcXD)s[\Vs\{XAXB)sJ = {[{siSi)Sc{s3S2)SD]s[\l\[{siS2)SA{s3S4)SB]sJ 

{Si S2 Sa 1 
S4 S3 Sb\. (20) 
Sc Sd S ) 

The matrix element of the operator f ^ = Sj ■ Sj can be derived similarly. For diagrams CI 
and C2, the matrix element is given by 

{iXcXD)s'J'^3\iXAXB)sJ = {[iSlS4)ScissS2)SD]s[h " | [(SiS2)S'A(s3S4)S's]fj 

= Sss'SsM-^f''^'''''''''''''SASBScSD 
Si S2 Sa 



X <! S4 S3 Sb\- [S.,{S,, + 1) - S,{S, + 1) - S,{S, + 1)] . 

The values of i, j, and Sij for diagrams CI and C2 are listed in Table I. 

Table I. The values of i, j, and Sij in Eq. (|2l|) . 



(21) 



Diagram 


i 


j 


Sij 


Prior CI 


1 


4 


Sc 


Prior C2 


2 


3 


Sd 


Post CI 


1 


2 


Sa 


Post C2 


4 


3 


Sb 



The matrix element of Vs = Sj-Sj for diagrams Tl and T2 is somewhat more complicated, 
and can be shown to be 
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{{XcXd)s[\Vs\{XaXb)sJ = {[{siSA)Sc{s3S2)SD]sMi ■ Sj\[{siS2)SA{s3SA)SB]s[) 
_^^5.4-.4-s.(25^3 + 1)(2^24 + 1)SaSbScSd 



Sis S24 
Si S2 Sa 

X < S3 S4 Sb 

Sl3 S24 S 




c 



s, 



S,iS, + l)-S,iS, + l)], (22) 



where i=l and j=3 for diagram Tl, and i=2 and j=4 for diagram T2. The quantities S13 
and 5*24 span the full allowed range in this summation. 

Eqs. (pOD, (pll), and (^) give the general results for the spin matrix element /spin of the 



unit operator and the Sj ■ s,- operator in a rearrangement collision. Our results agree with 



the explicit coefficients given in Table I of Barnes and Swanson 133 . 



VI. EVALUATION OF THE SPATIAL MATRIX ELEMENT 

In the quark-interchange reaction of Eq. the masses of the quarks and antiquarks 
are different in general. Previously, meson scattering calculations with unequal masses 
using this approach had been discussed in detail in coordinate space |3^. Here we give the 



corresponding momentum space results for general quark and antiquark masses. 
The spatial matrix element in Eq. ([TT|) is 

{{^C^D)fj'^-S'J^r\{^A^B)fj^-S.)) 
MaMbMqMu 

X {^c{LcMc)^D{LDMD)\Vr\^A{LAMA)^B{LBMB))). (23) 

For the four diagrams in the reaction A + B C + D , the spatial matrix element 

/space = {^c{LcMc)^D{LDMD)\Vr\^A{LAMA)^B{LBMB)) (24) 

can be written in the form 

= jj dK dK' ^AmkA - KA)]<^BmkB - Kb)] 
X $c[C(2fcc - Kc)]<l>DmkD - Kd)]V{k,' - k). (25) 



I space 



Here the momentum arguments are shown explicitly, and the angular momentum quantum 
numbers Lj and Mj for each meson are implicit. The quantity C = ±1 is an overall sign 
which depends on the diagram (see Table II). The quantity V{q), where q = k, — k', is the 
Fourier transform of the interaction Vijirij) [the spin and spatial part of Hijivij) in Eq. (|^)], 

V{q) = j dr e-*-^-^- \/,,(r,,). (26) 

The momenta n is the initial three-momenta of the scattered constituent that is initially in 
meson A and k' is its final three-momenta. The variables {fcj, {i = A, B, C, D)} are either 
K, or k' depending on the diagram, as specified in Table II. We shall use the bold-faced 
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symbols A, B, C, and D to represent the momenta of A, B, C, and D respectively. For 
simplicity we shall treat the scattering problem in the center-of-mass frame, so that B=—A 
and D=—C. The quantity {K^} is a function of A, C, and the mass parameter /j, which 
is a function of the quark and antiquark masses in meson i. The function Ki{A, C, fi) is 
tabulated for each diagram in Table II. For diagrams Tl and T2, the post and prior variables 
are identical and so do not need to be tabulated separately. 



Table 11. Diagram-dependent momentum arguments in post and prior formalisms. 



Diagram 


c 


kA 


Ka 




Kb 


kc 


Kc 


ko 


Kd 


Prior CI 


1 


K 


IaA 


n 


fBA + 2C 




fcC 




foC + 2A 


Prior C2 


-1 


K, 


IaA 


n 


fsA - 2C 


K 


-fcC + 2A 




-foC 


Post CI 


1 


K 


IaA 




fsA + 2C 


h! 


fcC 


k' 


foC + 2A 


Post C2 


-1 


K 


IaA 


h! 


fsA - 2C 


K 


-fcC + 2A 


K 


-foC 


Tl 


1 


K 


IaA 


k' 


fsA + 2C 


k' 


fcC 


K 


fnC + 2A 


T2 


-1 


K 


JaA 


k' 


IbA - 2C 


K 


-fcC + 2A 


k' 


-foC 



The mass parameter fi is unity for mesons with equal quark and antiquark masses. For 
unequal masses, the {fi} are tabulated in Table III in terms of 

^ ^ m{q)i - m{q)i 
m{q)i + m{q)i' 

The { fi} are the same in post and prior formalisms. 

Table III. The mass parameters fi for each diagram. 



Diagram= 


CI 


C2 


Tl 


T2 






1- A^ 




1- A^ 


fB = 


1- As 


l + As 


1- As 


1 + Ab 


fc = 


1 + Ac 


1 + Ac 


1 + Ac 


1 + Ac 


fD = 


1- Ad 


I -Ad 


1- Ad 


1- Ad 



To evaluate the spatial overlap integrals, we expand each meson wavefunction $(2p) as 
a linear combination of (nonorthogonal) Gaussian basis functions 0n(2p) of different widths 
as 



where 



N 



$(2p) = ^a„0„(2p), (28) 



n=l 



M^P) = ^n{2py^j^^ YUp) e^P{-|^}- (29) 

This expansion makes the spatial integrals tractable. We choose to normalize the basis 
function 0n(2p) according to 
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(ip|042p)r = l, (30) 
which leads to the normahzation constant 

_ / 1 y/' 1 

We also normalize the meson wavefunction $(2p) according to 

dp\^{2p)\^ = l, (32) 



which implies a constraint on the coefficients {a„}. 

We shall ffist present our results for the spatial matrix element Eq. (pSf ) for the case of 
all 5- wave mesons, each with a single Gaussian wavefunction of the type of Eq. (p^), 

(/.,[C(2/e. - Ki)] = N,exp{-j{2ki - K,)'} (33) 
where Aj = l/AnP'^ and 

^« = v^f7j ■ (34) 

The product of wavefunctions in Eq. (^) is explicitly 

0^[C(2fc^ - Ka)] (pBrn^B - Kb)] 0c[C(2fcc - Kc)] MCi^ko - Kd)] 

4 \ 

= NANBNcNnexp{-Y,^{2k, - K,f}. (35) 

1=1 ^ 

The argument of the exponential, from the product of the four meson wavefunctions, is a 
function of fcj = {k, k'} and the quantities {Ki}. It can also be written as a function of 
p = (k + k')/2 and q = k' — k. In terms of p and q, the {fcj} are given by 



where rji is 



ki=p-r]i q/2, (36) 



"■H!;; (37) 



Using Eq. (^Bj) and completing the square in the exponential, we obtain 

E J (2fc. - = 2 E A. (P - Po)^ + f (q - ^o)' + ^, (38) 

where the quantities po, Qq, fi, and u are defined below. The shift po is 

Po = roQ + So, (39) 
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where ro and Sq are 



^0 = h'^^i 



1=1 



1=1 



(40) 



(41) 



i=l 



i=l 



The quantity in Eq. (153) is 



(42) 



U=l 



the shift is 
2 

9o - 



\i=l 



u=l 



(43) 



and V is 



(44) 



i=l 



%=1 



The product of wavef unctions in Eq. (|25| ) can therefore be written in a shifted Gaussian 
form, 

0^[C(2fcA - Ka)] </)B[C(2feB - Kb)] 0c[C(2fcc - Kc)\ MCC^kn - Kd)] 



NaNbNcNo exp{-2 ^ A,(p - po)' " ^(^ " 5o)' - 



(45) 

The spatial matrix element of Eq. ( 125|) then becomes 

11 dK dK,' (j)AmkA - Ka)] 0B[C(2fcij - Kb)] 0c[C(2fcc - Kc)] MCC^kn - Kn)] V{q) 

4 

= lldpdqNANBNcNnexp{-2j2Kip-poY-^iq-qoY-'^}Viq). (46) 



The integration over p can be carried out analytically, which gives 

//rfp rfg n0aC(2fc. - K,)] V{q) = NaNbNcNd e'^ | rfq e^t^'^-'^o)^ \/(q). (47) 



Thus, the six-dimensional integral of Eq. ( ^61) is simplified to a three-dimensional integral 
involving an integration over V{q). 

The interaction V{q), which we take from the standard quark model Vij{r), is the sum 
of Fourier transforms of color- Coulomb, spin-spin contact, linear confinement and constant 
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terms. The sum of all diagrams arising from the constant term is zero for the scattering 
problem, so we do not need to include Vcon in deriving scattering amplitudes. 

For the remaining three interactions, we have (using some integrals of Ref. ||38|| ) 



(48) 



'dq e-^^'^-'^o? VM = [dq e't ('^"'^0)^ ( Z^] e-'^V4.^ 



3mimj\fiJ \l + 2cr^^ij "^""^^ 2{l + 2a^n)^' ^^^^ 
and 

/ dq e-tC'-o)^ VM = {-I) 8nb (f )''V^e--^/^F, (-i, | ^) , (50) 

where Vii^{q) = (—3/4) / dre~^'^'^br. These results allow one to evaluate the transition 
matrix elements Ty-j explicitly for the different interactions in diagrams CI, C2, Tl, and T2 
for the case of Gaussian meson wavef unctions. 

The wavefunctions we employ here are in general sums of Gaussians of different widths 
(Eq. (p8|)). Eq. (pSf ) can be evaluated in that case as well, so that the spatial matrix element 
Eq. ( p5D becomes 

-^space 

= jjdK dK,' <l>AmkA - Ka)] <fB[C(2fcB - Kb)] <fc[C(2fcc - Kc)] <fD[C(2fcD - Kd)] V{k' 

N N N N 

= X! X! X! X! (^nA(^nBO'ncO-no ^space{nA,nB,nc,nD), 
nyi=l n.B=l no=l n.£)=l 

where 

Ispace{nA,nB,nc,nD) = JJ dK, dn' 0^„^[C(2fc^ - Ka)] 0B„s[C(2fcB - Kb)] 

X 0cne[C(2fcc - Kc)] (l)DnM'^kD - K d)] V {k' - H,). (52) 

Ispace{nA,nB,nc,nc) is the previous result of Eq. (^) for a single component wavefunc- 
tion. The overlap integral in the multicomponent case is simply a sum of single-component 
contributions, each weighed by a coefficient product an^anj^anadnD- 

After the matrix elements for the interaction (^ are evaluated, the cross section for the 
process A + B — >■ C + D can then be obtained using conventional scattering theory, as 
discussed in Section II. 
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VII. EVALUATION OF THE SPATIAL MATRIX ELEMENT FOR AN X = 1 

MESON 



In the last section, we considered the scattering of S'-wave (L = 0) mesons. Here we 
generahze to colhsions in which a P-wave {L = 1) meson A colhdes with an S'-wave meson 
B, and scatter into two S'-wave mesons C and D. 

First we consider single-component Gaussian wavef unctions. (The results can be easily 
generalized to multicomponent Gaussian wavefunctions.) Eq. (p5D becomes 

II di^ dH,' M^kA - Ka) (psC^kB - Kb) (pcC^kc - Kc) (poC^kn - Kd) V{q) 

= lld,^ d^'NA \^Pa\''"^I^Z^^_'^LamAPa) exp{-^{2kA - Ka)'} 
X(PB{2kB - Kb) (pci2kc - Kc) (poC^kD - Kd) V{q), 
where 2pA = 2k a — Ka- Setting La = I for the P-wave meson A, we have 

I ^ r 2k, - JaA, if Ma= 

\2kA - KaI"-^ t/TTTr— "^l^mAPa) = -(2/^. - IaA,) - i{2Ky - fAAy) if Ma= 1 . (53) 
V ^2^^ + [ {2k, - JaA,) - ti2Ky - fAAy) if Ma=-1 



It then suffices to evaluate 



dn dK'{2K - f^A)kNAexp{-^i2kA - Ka)'} 

X <PB{2kB - Kb) (pcC^kc - Kc) 0D(2fcD - Kd) V{q) 

4 

dp dq {2k - JaA), NANBNcNDexp{-2Y,X,{p - Po)' - ^{q - q^)' - ^} V{q). (54) 

i=i ^ 

We can express 2k — A in terms of p and q; 

2k. - /aA = 2p-q- JaA = 2{p - po) + 2po-q- /aA . (55) 

Substituting Eq. (p9D into this result, we find 

2k - /aA = 2{p - Po) + (2ro - l)(q - q^) + (2ro - 1)^0 + 2so - /aA . (56) 

The integral of p — po gives zero. The integration over the last three terms can be carried 
out in the same way as in the = case, because Qq, Sq and A are independent of the 
integration variables. It is thus only necessary to evaluate the integral 

4 

dp dq NANBNcNDiq - q,)k exp{-2 Hp ' Pof Qo)' ' ^} V{q). (57) 

2 = 1 ^ ^ 



The integration over p can be carried out analytically, which gives 



4 



dp dq NANBNcNDiq - qo)fc exp{-2 ^ A,(p - po)^ - ^{q - q,)' - ^} V{q) 

i=i ^ ^ 

/ \ 3/2 a 

NaNbNcNd[—^A e-'i-—- f dq e^f ('^-'^o)V(q). (58) 
V2Ei=iA^; ^^^qokJ 
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The dq integrals for the various potentials have already been obtained in closed form, and 
the differentiation with respect to gofe is straightforward. We then find 



d 



Ok 



Jdq e-^^'^-'^of VcouXq) 



^^4' 2' 2 ^ 3 4' 2' 2 ' 



(59) 



d 



Ok 



J dq e-t('^-'^o)V,,(q) = 



— 87ra, / 27T 



3/2 



X 



V l + 2a2/i^ 



2aV ^ 
1 + 2(t2^^ 



3/2 



exp{- 



ml 



2(1 + 2a2/i) 



} 



(60) 



and 



d 



dqi 



Ok 



I dqe-i^"-kl'v,M 



.2).8.(f)'%v-^'= 

I 2' 2' 2 / 3 \2' 2' 2 



^qok- (61) 



The scattering amplitude Tji and cross section a/j for the dissociation of a P-wave meson 
through an SP — > SS" reaction will subsequently be evaluated using these results. 



VIII. MESON WAVE FUNCTIONS 

In nonrelativistic reaction theory, the equality of the scattering amplitude for rearrange- 
ment reactions in post and prior formalisms follows if and only if the two-body interaction 
used in evaluating the scattering matrix elements is identical to the interaction that gen- 
erates the bound state wavef unctions. If this is not the case, the post and prior scattering 
amplitudes will differ. It is therefore especially important to determine accurate bound state 
wavefunctions in evaluating scattering amplitudes. For this purpose, we will first search for 
a set of interaction Hamiltonian parameters that fit the known meson spectrum with reason- 
able accuracy. This interaction and the associated bound state wavefunctions will be used 
in our subsequent reaction calculations. 

For a bound state of quark i and antiquark j of momentum pi and pj and reduced 
momentum p, 

p= _\ (62) 

the Hamiltonian is 



rrii + rrij 



H=^ + Vir), (63) 

where is the reduced mass mimj/{mi + mj), and V{r) is the quark- antiquark interaction, 
(see i/,, of Eq. (|)), 
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For a quark and antiquark in a color-singlet hadron, the matrix element of —X{i) ■ X^{j)/4: 
is the familiar color factor Cf = —4/3. 

For given orbital angular momentum quantum numbers / and m, the eigenstate $(2p) of 
this Hamiltonian can be represented as the expansion in Eq. (^) in the set of (nonorthogo- 
nal) Gaussian basis states {0„} of Eq. (|29|) with expansion coefficients {an}. The eigenvalue 
equation iJ$ = then becomes the matrix equation 

Ha = EBa, (65) 

where a is a column matrix with elements {ai, 02, oat}, B is the symmetric matrix 

/ 1+3/2 

B,, = ^Jdp cl)*{2p) 0,(2p) = ( j , (66) 
and Ti. is the Hamiltonian matrix 

n^J = Tij + (67) 
which is the sum of the kinetic energy matrix Tij 

T,, = /* A'(2P) I ^,(2p) = (21 + 3)-f- (-^j 1^ . (68) 

and the potential matrix Vij 

V,, ^ {t\V\j) = i2nfl dr 0*(r) \/(r) 0,(r), (69) 
where 0j(r') is the Fourier transform of 0j(2p), 

(^^"■>^(2P)- (70) 
Evaluation of the potential matrix elements for the color-Coulomb interaction gives 
VijiCou.) = Cf4:7ias{i\ ^ \j) 

= ^^ (2vr)3/M2/ + l)!! V^ + H^J ' ^ ^ 

For the linear interaction we similarly find 

y,,(lin.) = (z|C;(-^)6r |j) 

1+3/2 



r 8^ (^ + 1)' 1 Cblll\ (72) 

47 /5 (27r)3/2(2/ + l)!! v/rrj l^z + jj ' ^ ^ 
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and for the spin-spin interaction we find 

,-3 



V^,(ss) = -C 



a 



Smirrij 



IT' 



3/2 



Z+3/2 



Given these Hamiltonian matrix elements, the eigenvalues and eigenvectors can be obtained 
from the eigenvalue equation (^). In our numerical calculations we used a six-component 
[N = 6) space of Gaussian basis functions. 
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Fig. 3. Comparison of experimental energy levels and 
theoretical energy levels calculated with the parameter 
set of Eq. (|74D. 

For this study we assumed a running coupling constant combined with an otherwise 
fairly conventional quark model parameter set, given by 

12n 



asiQ 



A =10, B = 0.31 GeV, 



(33 - 2n;) ln(A + ■ 
6 = 0.18GeV^ or = 0.897 GeV, m„ = = 0.334 GeV, 
m, = 0.575 GeV, = 1.776 GeV, mb = 5.102 GeV, 
Vcon = -0.620 GeV. (74) 

We identified the scale in the running coupling constant with the square of the invariant 
mass of the interacting constituents, = Sij. This set of parameters leads to a meson 
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spectrum which is reasonably close to experiment (see Fig. 3); the theoretical masses and 
wavefunctions are given in Appendix A. The parameter set used earlier in ^| is similar 
to this set but it has a fixed strong coupling constant. An alternative set of quark model 
parameters, without a running coupling constant, was used for comparison. This second 
set was as = 0.594, b = 0.162 GeV^ a = 0.897 GeV, m„ = = 0.335 GeV, and rric = 
1.6 GeV, and a fiavor- dependent Vcon- 

Having obtained the set of coefficients {an} for each initial and final meson, we can 
proceed to the calculation of the scattering amplitudes Tjj and the dissociation cross sections 
afi. Some explanation of the evaluation of our (somewhat arbitrary) choice of momentum 
scale = Sij in the running coupling constant <y{Q'^) in Eq. (0) for the scattering problem 
is appropriate. For a reaction process involving the interaction of constituent i = a in meson 
A{aa') with j = b in meson B{bb'), we can determine the invariant mass squared Sij of a and 
b as follows: Constituent a carries a fraction x+ of the forward light-cone momentum of the 
initial meson A, and b carries a fraction of the backward light-cone momentum of initial 
meson B. For simplicity, we assume that the momentum fraction carried by a constituent 
is proportional to its rest mass, which is exact in the weak binding limit; 



x+ 



rria + nia' 



(75) 



X- 



rrib 



nib + my 

Assuming also that constituents a and b are on mass shell, their momenta are 



(76) 



and 




x+{Aq + A,) ± 



x(B.-B^± 



(Ao + A, 



mt 



{Bo — Bz 



(77) 



(78) 



(79) 



and the invariant mass of a and b is then given by 

Sab = {a + bf = (ao + bof - {a^ + b^f. 

We identify this with the argument of the running coupling constant asiQ"^) in Eq. ([f4|). 

In Fig. 4 we show a test of the accuracy of this scattering model with experimental 
data in an analogous low-energy reaction, J = 2 vrvr scattering. The prediction for the 
dominant S'-wave phase shift is shown together with the data of Hoogland et al. [^. Note 
that this is not a fit; all the parameters are determined by meson spectroscopy, which fixes 
the interquark forces and wavefunctions that are then used to calculate the meson-meson 
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scattering amplitude. 
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Fig. 4. Comparison of our theoretical / = 2 tttt phase 



shift (solid curve) with the data of Hoogland et al. [41 



IX. A TEST OF 'POST-PRIOR' EQUIVALENCE IN vr+J/V^ DISSOCIATION 

Fig. 1 shows the 'prior' meson-meson scattering diagrams, in which the interactions occur 
before constituent interchange. There is a corresponding 'post' set of diagrams, in which 
the interactions occur after constituent interchange, as shown in Fig. 2. As we noted earlier, 
these are equivalent descriptions of the transition matrix element, and should give identical 
cross sections. 

In nonrelativistic scattering theory the post and prior results can be formally proven to 
be equivalent, so that the two theoretical cross sections are indeed identical. This proof 
requires that the interaction V used to determine the external meson wavefunctions be 
identical to the interaction used in the evaluation of the scattering amplitude. Numerical 
confirmation of this post-prior equivalence constitutes a very nontrivial test of the accuracy 
of the numerical procedures used both in determining the bound state wavefunction and 
in evaluating the complete meson-meson scattering amplitude. This post-prior equivalence 
was discussed in detail and demonstrated numerically by Swanson for tttt pp scattering 
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34|, where its relevance to establishing Hermitian effective scattering potentials was 



in Ref. 
shown. 

To test post-prior equivalence in our J / ip dissociation calculations (in the nonrelativistic 
formalism), we have carried out the evaluation of the cross section using both post and prior 
formalisms. Of necessity we assumed nonrelativistic kinematics and theoretical masses to 
determine the external meson momenta |A| and |C|, which appear in the overlap integrals. 
The post-prior equivalence holds if the coupling constant does not depend on energy. We 
are well advised to use a set of parameters with a fixed running coupling constant for 
this test. We used therefore the parameter set |2T| = 0.58, 5 = 0.18 GeV^, (j=0.897 GeV, 
^u,d = 0.345 GeV, nic = 1.931 GeV, and Kon = —0.612 GeV, which are close to standard 
values and were chosen because they reproduce the physical masses of the vr, D, D*, J /ip and 
rather well. Fig. 5a shows the dissociation cross sections for ir+ip collisions as a function 
of the initial kinetic energy E^e measured in the center-of-mass frame, defined as E^e = 
\fs — TTiA — rnB where rriA and ttib are the masses of the incident mesons. The differences 
between the post and prior results in Fig. 5a are indeed rather small, which confirms post- 
prior equivalence assuming non-relativistic dynamics. (The residual discrepancy is mainly 
due to our use of a truncated basis in expanding the meson wavef unctions.) Fig. 56 shows 
the corresponding results for n+ip' dissociation. In this case there is greater sensitivity to 
the approximate ip' wavefunction, due to large cancellations in the numerical integration of 
a radially-excited wavefunction. 
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Fig. 5. Comparison oiTT + J/ip and Tc + ip' cross sections 
derived using post and prior formalisms. 



In the comparisons shown in Fig. 5 we have used theoretical masses for the mesons; 
the proof of post-prior equivalence makes use of the theoretical bound-state masses from 
the Schrodinger equation with the given potential, rather than the experimental ones. Our 
theoretical masses are close to experiment but are not exact, as is evident in Fig. 3. If 
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one instead assumes experimental values for the meson masses, there will be additional 
post-prior discrepancies in our cross section calculations. In our subsequent cross section 
calculations we do assume experimental masses in order to reproduce correct thresholds; the 
consequence is a systematic uncertainty in the cross sections which may be estimated by 
comparing post and prior predictions. 

It should be noted that 'post-prior equivalence' had only been proven for bound-state 
scattering in non-relativistic quantum mechanics Recently an extension of this proof 



to scattering in a relativistic generalization of quantum mechanics was given by Wong and 
Crater using Dirac's constraint dynamics. A full relativistic treatment of this problem 
will involve the derivation of relativistic two-body wavefunctions and Wigner spin rotations, 
which is beyond the scope of the present investigation. We will adopt an intermediate ap- 
proach and assume relativistic kinematics for the initial and final mesons and use relativistic 
phase space; in consequence we find different 'post' and 'prior' cross sections in general. Here 
we will take the mean value of the 'post' and 'prior' results as our estimated cross section, 
separate 'post' and 'prior' cross sections will be shown as an indication of our systematic 
uncertainty. 



X. CROSS SECTIONS FOR J/V' AND V' DISSOCIATION 

Depending on the incident energies, dissociation of the J /if) and ip' by pions can lead 
to many different exclusive final states. There are several selection rules that eliminate 
or suppress many of the a priori possible final channels, given our simple quark-model 
Hamiltonian and Born-order scattering amplitudes. Considerable simplification follows from 
the fact that our model Hamiltonian conserves total spin Stot- Since the J /if) and ip' have 
S = 1 and pions have S" = 0, the initial and final states in t^+J /ip and tt+iIj' collisions must 
both have Stot = 1; this forbids DD final states. With increasing invariant mass we next 
encounter the final states DD*, D*D and D*D* . In Fig. 6 we show the exclusive dissociation 
cross sections for -K+J/ip and ir+ip' collisions into these first few allowed final states. The 
total dissociation cross section, which is the sum of the exclusive cross sections, is shown as 
a solid line. Our estimate is the mean of the total 'post' and 'prior' cross sections, which are 
also shown in Fig. 6. The estimated range of uncertainty, due to the post-prior discrepancy 
and to parameter variations, is shown as a shaded band. 

The J/ip dissociation processes ti+J/iIj D*D or DD* have a threshold of about 
i?x_B=0.65 GeV and the total dissociation cross section reaches approximately 1 mb not 
far above threshold (Fig. 6a). This tt + J/ip cross section is considerably smaller than the 
peak value of about 7 mb found by Martins et al., largely due to their assumption of a 
color-independent confining interaction. 

The threshold for it + tp' DD* + D*D is only about Eke = 0.05 GeV, and in con- 
sequence the ip' cross sections are much larger near threshold. The total ir+ip' dissociation 
cross section has maxima of ~ 4.2(0.5) mb and ~ 2.8(0.5) mb at Eke ~ 0.1 GeV and ~ 0.22 
GeV respectively (Fig. 66). It is interesting that the exclusive ip' dissociation cross sections 
are very small near Eke = 0.3 (for the DD* final state) and 0.4 GeV (for the D*D* final 
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state), due to almost complete destructive interference between the diagrams. 
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Fig. 6. Exclusive and total dissociation cross sections for 7r+ J /ijj 
(Fig. 6a) and Tv+ip' (Fig. 66). In each panel the solid curve gives 
our estimated total cross section, which is the mean of the 'prior' 
and 'post' results. 



The relative importance of the various terms in the Hamiltonian in these dissociation 
amplitudes is of course a very interesting question. Unfortunately it is also somewhat am- 
biguous, because the individual terms differ between post and prior formalisms; only the 
sum of all terms is formalism independent in non-relativistic quantum mechanics. We find 
that the spin-spin interaction makes the dominant contribution to n+J/ip dissociation in 
the prior formalism; ir+ip' dissociation in contrast is dominated by the linear confining in- 
teraction. In the post formalism we find that both Tc+J/ip and n+tp' are dominated by the 
spin-spin interaction. In all these cases the color-Coulomb contribution is subdominant. 

Our results have interesting consequences for the survival of J/ip and ip' mesons propa- 
gating in a gas of pions. The pions produced in a heavy ion collision have a roughly thermal 
distribution, with T ^ 200 MeV in the nucleus-nucleus center-of-mass system, whereas heavy 
quarkonia such as the J/ip and ip' are produced approximately at rest. Thus the relative 
kinetic energy is typically a few hundred MeV. This is below the n+J/ip dissociation thresh- 
old, but above that of n+ip', and in consequence we expect n+ip' collisions to deplete the ijj' 
population much more effectively than n+J/ip collisions reduce the initial J/ip population. 
The weakness of n+J/ip dissociation is further reduced by the small cross section we find 
for the n+J/ip relative to n+ip'. 

Next we consider the very interesting exothermic collisions of charmonia with light vector 
mesons, specifically p+J/tp and p+ip'. Since the p meson has S = 1, the total spin of the 
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p + {J l'>p or ip'} system can be Stot = {0, 1, and 2}. This Stot is conserved in our model, 
and so must agree with the Stot of the final state. The low-lying open charm final states are 
DD with Stot = 0, DD* and D*D with Stot = 1, and D*D* with Stot = {0, 1, and 2}. The 
contribution of transitions to radially- and orbitally-excited final states will be considered 



in future work 43 
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Fig. 7. The total unpolarized dissociation cross sections are 
shown as the solid curves for p+J/ijj in Fig. 7a and for p+ip' in Fig. 
7b. Unpolarized exclusive dissociation cross sections a'j^'^^^lStot) for 
different final states and different 5*40* are also shown. 



In the collision of an unpolarized p with an unpolarized J/ip, the total dissociation cross 
section comprises of contributions {Stot) from different final states / and different total 
spin values Stot of the system, 

<Ttr' = T.T.^r"'(Stot), (80) 

/ Stot 

where for this case with La = and 5^ = 7^ we can deduce from Eq. (|T7|) the following 
relationship 
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(81) 



and crf{Stot) is the dissociation cross section for the final state / when the initial two- 
meson system is prepared with a total spin Stot- [In our earlier work for p -\- J/ijj 
and p + ip' collisions, <7f{Stot) results were presented and the total cross section of atot = 
J2f Y^Stot ^i^tot) was evaluated. However, for the collision of unpolarized mesons, one should 
use the unpolarized total dissociation cross section as given by Eqs. ( pOD and (|8T|).]. 

The unpolarized total p + J/ip dissociation cross section are shown in Fig. 7a. The 
exclusive cross sections (TY''^"^{Stot) for dissociating into various final states in an unpolarized 
collision are also exhibited. The reaction p + J/iIj ^ DD is exothermic, so the cross section 
(J^£)°^{Stot = 0) diverges as l/y/ExE as we approach threshold. For other p + J/ip exclusive 
final states the thresholds lie at sufficiently higher energies to be endothermic. We find 
an unpolarized total p + J/ip dissociation cross section of 2.4(0.4) mb at Eke = 0.1 GeV, 
which has decreased to about 1.9 mb at Eke = 0.2 GeV. At very low kinetic energies, the 
contribution to the dissociation of J/ip hj p comes from the DD{Stot = 0) final state. At 
slightly higher energy, it comes mainly from DD* and D*D{Stot = 1) final states. At Eke 
about 0.2 GeV, it comes dominantly from the D*D*{Stot = 2) final state. 

We have carried out similar calculations for p + ip' collisions, and the results are shown in 
Fig. 7b. In this case the reactions p + ip' ^ DD, DD*, D*D and D*D* are all exothermic, so 
all these exclusive cross sections o"j"^°^(S'tot) diverge as I/^/Eke as we approach threshold. 
The dominant contribution to the dissociation comes from the D*D*{Stot = 2) final state. 
The total unpolarized p + ip' dissociation cross section falls from 4.5(0.1) mb at Eke = 
0.1 GeV to 1.5(0.3) mb at Eke = 0.2 GeV and 0.4(0.2) mb at Eke = 0.3 GeV. 

In p+J/ip dissociation the dominant scattering contribution in the prior formalism is 
due to the linear interaction. In the post formalism the dominant contribution arises from 
the color-Coulomb and linear interactions at energies Eke < 0.1 GeV, from color-Coulomb 
at 0.1 GeV< Eke <0.4 GeV, and from the spin-spin interaction at Eke > 0.4 GeV. 

We turn next to the dissociation of J/ip and ip' in collision with K. Our predictions 
for K + J/ip and K + ip' dissociation cross sections are shown in Fig. 8. The K + J/ip 
process has a threshold kinetic energy of about 0.4 GeV, and the maximum cross section is 
about 0.7 mb. In K+ip' dissociation the reactions are exothermic for the allowed final states 
Ds + D* , D* + D and D* + D*. The total K+ip' dissociation cross section, shown in Fig. 86, 
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is about 1 mb at Eke ~ 0.4 GeV, and also diverges as 1/\/Eke as we approach threshold. 
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Fig. 8. Dissociation cross sections for K+J/ip (Fig. 8a) and K+ip' 
(Fig. 86). 

The model of Barnes and Swanson has thus far received experimental support from 
extensive comparisons with light hadron scattering data {1 = 2 nn ||3^, 7 = 3/2 Kn ||35|| , 
/ = 0, 1 KN 1^61 )• '^^^ results of J/tp and ijj' dissociation in collision with vr, p, and 



K need to confront experimental data through an examination of their implications on J/ ip 
and ip' suppression in high-energy heavy-ion collisions. The extent to which the observed 
anomalous suppression in Pb+Pb collisions is due to the dissociation J / ip by vr, p, and K will 
require further quantitative study. J/ip suppression by these mesons must be incorporated 
and subtracted in order to identify the suppression of J/ip production by the quark-gluon 
plasma. It is therefore important to carry out a detailed simulation of J/ absorption using 
cross sections obtained here in order to understand the nature of J/ip suppression in Pb+Pb 
collisions and to provide a test of the theoretical J/ ip dissociation cross sections. 

XI. CROSS SECTIONS FOR T AND T' DISSOCIATION. 



It has been noted that the suppression of the production rates of bb mesons, such as the T 



and T', may also be useful as a signal for QGP production (see Ref. and references cited 
therein). It is of interest to calculate the T and T' dissociation cross section in collisions 
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with TT, p, and K so as to infer the effects of T and T' suppression by hadronic produced 
particles. 
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Fig. 9. The dissociation cross sections for T (Fig. 9a) and T' (Fig. 
96) in colhsion with vr. Note that the scales of Figs. 9a and 96 are 
different . 



We show in Fig. 9a the total cross section for the dissociation of the T in collision with 
71. The threshold is at Eke ~ 1 GeV, and the maximum cross section is about 0.05 mb. 
The small cross section arises from the combined effects of a large threshold and a small 
value of the strong interaction coupling constant. We show the dissociation cross section for 
T' in collision with vr in Fig. 96. It has a threshold of E^e ~ 0.45 GeV, and the peak cross 
section is about 5 mb at Eke ~ 0.55 GeV. This dissociation cross section is much larger 
than that for the dissociation of T by vr. 

In Fig. 10a we show the unpolarized cross section for the dissociation of T in collision 
with p. Unpolarized exclusive dissociation cross sections (J^^"^"^ (Stot) for different final states 
and different Stot are also shown. The reaction process is endothermic with a threshold at 
Eke ~ 0.3 GeV and a peak cross section of 0.15 mb at Eke ~ 0.45 GeV. In Fig. 106 we 
show the dissociation cross section for T' in collision with p. The reactions are exothermic 
and the total dissociation cross section behaves as I/^/Eke near Eke ~ and decreases to 
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about 1.6 mb at Eke ~ 0.2 GeV. 
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Fig. 10. Unpolarized total dissociation cross sections and unpolar- 
ized exclusive dissociation cross sections cr"°^°'(5'tot) for p + T (Fig. 
10a) and p + T' (Fig. 106) for various channels and total spin Stot- 
Note that Figs. 10a and 106 have different scales. 



In Fig. 11a we show the cross section for the dissociation of T by K. The threshold is at 
Eke ~ 0.75 GeV, with a peak total cross section of about 0.05 mb at Eke ~ 0.85 GeV. We 
show in Fig. 116 the dissociation cross section of T' by K. The threshold is at Eke ~ 0.2 
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GeV, with a peak total cross section of about 2 mb at Eke ~ 0.25 GeV. 
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Fig. 11. The dissociation cross sections for K + T' (Fig. 11a) and 
K + T' (Fig. 116) for various final states. Note that Figs. 11a and 
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XII. CROSS SECTIONS FOR XJ DISSOCIATION 

We can calculate the dissociation cross sections of xj mesons in collision with vr, p, and 
K using the quark- interchange model. A xj meson has a spin quantum number S = 1 and a 
TT has S = 0. The collision of a xj with a vr gives rise to a system with a total spin Stot = 1- 
On the other hand, the interaction of Eq. (2), which leads to the dissociation, conserves the 
total spin. Therefore the lowest-energy final states are DD*, D*D, and D*D*. 

We show the results for the dissociation cross sections of unpolarized xj m collision with 
TT in Fig. 12. The unpolarized dissociation cross sections for the final states of DD* and 
D*D* are shown as the dotted and the dashed-dot curves respectively. The unpolarized total 
cross sections for scattering into these lowest channels are shown as the solid curve obtained 
as the mean value of the total cross section from the 'prior' and 'post' formalisms. The 
estimated systematic uncertainty in the total cross section due to the post-prior discrepancy 
is again indicated as a band in the figure. The dissociation of xo by vr has a threshold 
of Eke ~ 0.32 GeV and a peak cross section of 1.53 mb at Eke ~ 0.5 GeV (Fig. 12a). 
The dissociation of Xi by tt has a threshold of Eke ~ 0.23 GeV and a peak dissociation 
cross section of 2.41 mb at Eke ~ 0.46 GeV (Fig. 126). The dissociation of X2 by vr has 
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a threshold of Eke ~ 0.18 GeV and a peak dissociation cross section of about 2.98 mb at 
Eke ~ 0.41 GeV (Fig. 12c). 
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Fig. 12. Unpolarized total cross sections and exclusive unpolarized 
cross sections for the dissociation of xj into D + D* and D* + D* 
in collision with vr. Fig. 12a is for it + Xo, Fig- 126 for vr + xi, and 
Fig. 12c for vr + X2- 



It is interesting to note that the maximum of the unpolarized total dissociation cross 
section for tt + X2 is only slightly greater than that for n + xi but is nearly twice as great as 
the maximum of the dissociation cross section for tt + Xo- This indicates that the dissociation 
of Xj is very sensitive to the threshold value. We found numerically that if the threshold 
value for n + xo were taken to be the same as the threshold value for ti + X2, the unpolarized 
dissociation cross sections would be the same. 

The dissociation amplitudes of the XjJz mesons in collision with pions depend on J^. A 
detailed discussion of the dissociation cross section for various J and Jz substates will be 



presented in |H3[. The dependence on is however quite weak. For the same value of J, 



dissociation cross sections of xjj^ in collision with n vary only by a few percent for different 
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The thresholds for vr + xj dissociation he in between those of ir + J/ip and tt + ip', and 
the maxima of the total dissociation cross sections for the n + Xj collisions are greater than 
that for the t: + J/ip collision but less than for the 7i + ip' collision. 




Fig. 13. Unpolarized total cross sections and unpolarized exclusive 
cross sections <j'^^'^°^ (Stot) for the dissociation of xj into DD, D+D* 
and D* + D* in collision with p. Fig. 13a is for p + Xoy Fig- 13^ for 
p + Xi, and Fig. 13c for p + X2- 



We show in Fig. 13 the dissociation cross sections of unpolarized xj in collision with 
unpolarized p. The lowest-energy final states are DD, DD*, D*D, and D*D* , characterized 
by different values of the total spin Stot- The unpolarized total dissociation cross section 
comprises of contribution o"j°^°'(S'tot) from different final states {/} and different total spins 
{Stot} of the system, 

cTtT' = Y.Y.^r''\Stot\ (82) 

/ Stot 

where for this case with = 1 and Sb = 1, crTP°^(S'tot) is more complicated than the 
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expression of Eq. (^) for La = and 5*^ = 1. It can be determined from Eq. (|T^) and is 
given by 

(Sa Sb S 

ar°'(^tot) = E i^jy {La La\ a{LAMASS,), (83) 
J^'^ [Ja Sb J ] 

where ^(L^M^S'S'^) is the cross section for the initial meson system to have a total internal 
orbital angular momentum La with azimuthal components Ma and total spin S. We show 
in Fig. 13 the unpolarized total dissociation cross section for p + xo in Fig- 13a, p + Xi in 
Fig. 136, and p + X2 in Fig. 13c. The exclusive cross section crj°^°'(S'tot) for different final 
states and Stot are also shown. The dissociation of xj in collision with p is exothermic. The 
dissociation cross sections have the common features that they diverge as 1/ \/Eke near 
Eke ~ and decreases monotonically as Eke increases. The dominant contribution to 
the dissociation cross sections comes from the D*D*[Stot = 2) final state. The unpolarized 
total dissociation cross section for p + xq at Eke = 0.05, 0.1, and 0.15 GeV are 8.0 mb, 3.5 
mb, and 1.6 mb respectively (Fig. 13a). The unpolarized total dissociation cross section for 
p + Xi at Eke = 0.05, 0.1, and 0.15 GeV are 5.5 mb, 2.0 mb, and 0.8 mb respectively (Fig. 
136). The unpolarized total dissociation cross section for p + X2 at Eke = 0.05, 0.1, and 
0.15 GeV are 4.3 mb, 1.5 mb, and 0.5 mb respectively (Fig. 13c). Thus, for the same kinetic 

energy Eke, (rtot"'\p + Xo) > o-to7°^(p + Xi) > (^tot"'\p + X2)- 

We show in Fig. 14 the unpolarized dissociation cross section of xj in collision with 
K. The lowest-energy final states are DgD*, D*D, and D*D*. For the dissociation of Xo 
in collision with K, the reaction has a threshold at 0.07 GeV. The total dissociation cross 
section rises from the threshold to a maximum cross section of 1.7 mb at 0.27 GeV (Fig. 
14a). For the dissociation of xi and X2 in collision with K, the reactions K + xi and K + X2 
are exothermic for the final states of DgD* and D*D. The total dissociation cross sections 
of Xi and X2 in collision with K behave as 1/\/Eke near Eke ~ but decrease very rapidly 
as Eke increases. These cross sections then maintain a relatively constant value up to 
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Eke ~ 0.2 GeV before decreasing at higher kinetic energies. 
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Fig. 14. Unpolarized cross sections for the dissociation of xj 
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XIII. DISCUSSION AND CONCLUSIONS 

We have used the Barnes- Swanson quark-interchange model, with parameters taken from 
fits to meson spectroscopy, to evaluate the low-energy cross sections for the dissociation of 
the J/V', ip' 1 Xi ^) s-iid T' in collision with vr, p and K. The cross sections obtained here 
should be useful as estimates of the importance of "comover" scattering in suppressing heavy- 
quarkonium production, which is of considerable interest in the search for the quark-gluon 
plasma. 

Our results show that the cross section for the dissociation of J/ip hy ti occurs at a 
relatively high threshold, and the peak total cross section is about 1 mb. In contrast, the 
cross section for the dissociation of ip' by tt occurs at a low threshold and the cross section is 
much larger. We have also evaluated the corresponding cross section for the dissociation by 
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p. This process is exothermic, and has a large total dissociation cross section that diverges 
at threshold. 

We previously noted that our n+J/ip cross section is considerably smaller than the 
estimate of Ref. ||2^, although we use a similar approach. There are several differences 
between the two approaches which lead to this discrepancy. First, Martins et al. assumed 
that the confining interaction is an attractive Gaussian potential which acts only between 
quark- ant iquark pairs. The neglect of the quark-quark and antiquark-antiquark confining 
interaction amounts to discarding the transfer diagrams (Tl and T2) for the confining po- 
tential. Since we find that the transfer and capture diagram confinement contributions are 
similar in magnitude but opposite in sign (due to color matrix elements), the confining inter- 
action assumed by Martins et al. did not incorporate this important destructive interference. 
Second, their use of a Gaussian, rather than the usual linear confining potential, obviously 
leads to quantitatively different cross sections. 

The destructive interference between transfer and capture diagrams with spin- 
independent forces (color- Coulomb and confinement) has been noted previously. (See, for 
example, Refs. p3| , p^ and references cited in [0.) This interference explains the well- 
known spin-spin hyperfine dominance in light hadron scattering in channels such as / = 2 
TTTT, and the core NN interaction. With heavy quarks, however, the hyperfine interaction 
contribution is smaller due to the large charm quark mass; for this reason we included 
the color-Coulomb and confining interactions in our analysis. Our results indicates that 
the spin-spin, the linear interaction, and the color-Coulomb interactions all give important 
contributions to the dissociation cross sections. 

It is of interest to compare our results of the dissociation cross section with those obtained 



in the meson exchange model with effective Lagrangians ||2^- p7|] . In the effective Lagrangian 
approach, the dissociation cross section increases with energy, as expected for the t-channel 
exchange of a spin-one particle. For example in Ref. pGf the dissociation cross section is 
about 30 mb for tt + J/^Ij ^ D + D* and about 80 mb for vr + T ^ 5 + B*, at 1 GeV 
above the threshold. These large cross sections continue to increase with increasing energy. 
In contrast, in our quark models calculation using the Barnes-Swanson model, we find very 
small cross sections this far above threshold for n + J/ip ^ D -\- D* and it + T ^ B + B*. 
The predicated peak cross section for 7i -\- J/tp ^ D + D* is about 0.5 mb and occurs at 
about 0.05 GeV from the threshold. The predicted peak cross section for vr + T — > 5 + S* 
is even smaller (about 0.03 mb), and it occurs at about 0.02 GeV from the threshold. These 
cross sections decrease rapidly at higher energies. Thus, the large cross sections obtained 
in the effective Lagrangian approach differ by orders of magnitude from the quark model 
results obtained here. We believe that the large increase in these dissociation cross sections 
predicted by the effective Lagrangian meson exchange models at high energies is unrealistic, 
since the momentum distributions of the boosted final and the initial states have little 
overlap at high energies. 

As the effective Lagrangian approach does not contain information about the internal 
structure of the interacting hadrons, phenomenological form factors have been introduced 
to reduce the theoretical cross sections [^5|-p7|. A realistic description of the form factors 
should incorporate the meson wave functions and the dynamics of the scattering process. 
Without a derivation of these form factors, one encounters considerable uncertainty, as 
experimental data on these reaction processes are unavailable. The form factors introduced 
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in 



25|-p7[] lead to changes of the theoretical cross section at high energies by several orders 
of magnitude. The results are sensitive both to the assumed coupling strength and to 
the functional dependence of the form factor. In view of the strong dependence of the 
theoretical results on the form factor and the coupling constants, a careful determination of 
these quantities are required in future work. 

Although there is no direct experimental measurement of these cross sections to which 
we can compare our results, the small vr+J/?/' and the large 7r+ip' dissociation cross section 
at low kinetic energies obtained here appear qualitatively consistent with earlier results in a 
microscopic model of J/V' and ip' suppression in 0+A and S+U collisions Hopefully, 
future Monte Carlo simulations of the dynamics of charmonium in heavy-ion collisions will 
lead to a more direct comparison. It is interesting to note that dissociation of J /iphj and p 
populate different states (for example, -n+J /ip does not lead to DD in our model but p+J/ip 
does). It may provide a way to separate these processes by studying the relative production 
of DD, D*D, DD*, and D*D*, if the open charm background can be subtracted. However, 
the large ratio of initial open charm to J/ip production in a nucleon-nucleon collision may 
make this separation very difficult. 

In the future it may be useful to carry out detailed simulations of J/tp absorption in 
heavy-ion collisions using the cross sections obtained here. If our cross sections do prove to 
be reasonably accurate, it will clearly be useful to incorporate them in simulations of hadron 
processes in relativistic heavy-ion collisions that use the J/ ip suppression as a signature of the 
quark-gluon plasma, in order to isolate the effects of J/ip suppression due to its interaction 
with hadron matter. 
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Appendix A: Tabulation of Bound State Wave Functions 

The wavefunction in reduced momentum 2p = pg — pg is represented as a linear combi- 
nation of Gaussian wave functions with linearly-spaced jS'^, of the form 

where $(2p) and 0„(2p) are normalized according to Eqs. (|3^) and (^OD, 

J dp m2p)\' = jdp |0„(2p)p = L (84) 

The coefficients {«„} for each meson in an A^=6 basis, with a different (3 for each meson, 
are listed in the following table. 
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